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Abstract 

We study two equations of Lotka-Volterra type that describe the Darwinian evolution of a 
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■ population density. In the first model a Laplace term represents the mutations. In the second one 

(-H we model the mutations by an integral kernel. In both cases, we use a nonlinear birth-death term 

that corresponds to the competition between the traits leading to selection. 

In the limit of rare or small mutations, we prove that the solution converges to a sum of moving 
Dirac masses. This limit is described by a constrained Hamilton- Jacobi equation. This was already 
proved in [8j for the case with a Laplace term. Here we generalize the assumptions on the initial 
data and prove the same result for the integro-differential equation. 
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O ■ 1 Introduction 



We continue the study, initiated in [8] , of the asymptotic behavior of Lotka-Volterra parabolic equa- 
tions. The model we use describes the dynamics of a population density. Individuals respond differently 
to the environment, i.e. they have different abilities to use the environment resources. To take this 
fact into account, population models can be structured by a parameter, representing a physiological 
trait, denoted by x € M. d below. We denote by n(t, x) the density of trait x. The mathematical 
modeling in accordance with Darwin's theory consists of natural selection and mutations between the 
traits (see [181 12H I2Z1 [25] for literature in adaptive evolution) . We represent the growth and death 
rates of the phenotypical traits by R(x,I). The term I(t) is an ecological parameter that corresponds 
to a nutrient that itself depends on the population. We use two different models for mutations. The 
simpler way is to represent them by a Laplacian. 

(d t n e -eAn e = ^R(x,I e (t)), x G M d , t > 0, 
1n e (t = 0)=n°GL 1 (M d ), n° e > 0, 
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I e (t) = / ip(x)n e (t,x)dx. (2) 

JR d 

Here e is a small term that we introduce to consider only rare mutations. It is also used to re-scale 
time to consider a much larger time than a generation scale. 

A more natural way to model mutations is to use an integral term instead of a Laplacian. 

(d t n e = ^R(xJ e (t)) + ±J±K(M^)b(y,I e )n € (t,y)dy, x G M d , t > 0, 
\n e (t = 0) = n° G L 1 (M d ), n° > 0, 

7 e (t) = / n e (t,x)dx. (4) 

Such models can be derived from individual based stochastic processes in the limit of large popula- 
tions (see plfTi]). 

In this paper, we study the asymptotic behavior of equations ([H)-© and ([3])-([2]) when e vanishes. 
Our purpose is to show that under some assumptions on R(x,I), n e (t,x) concentrates as a sum of 
Dirac masses that are traveling. In biological terms, at every moment one or several dominant traits 
coexist while other traits disappear. The dominant traits change in time due to the presence of mu- 
tations. 

We use the same assumptions as [8]. We assume that there exist two constants ip m , ipM such that 

< < < V>M < oo, V G W 2 ^{R d ). (5) 
We also assume that there are two constants < I m < Im < oo such that 

min R(x,I m ) = 0, max R(x, Im) = 0, (6) 

and there exists constants Ki > such that, for any x G I G R, 

B /? 

-K 1 < w (x,I)<-K^ 1 <0, (7) 

sup || R(-,I) || W /2,oo (R d)< K 2 . (8) 
We also make the following assumptions on the initial data 

I m < I ip(x)n° £ (x) < I M , and 3 A, B > 0, n° € < e~ A ™ +B . (9) 

Here we take ip(x) = 1 for equations ([3])-([4]) because replacing n by tpn leaves the model unchanged. 
For equation (|3|) we assume additionally that the probability kernel K(z) and the mutation birth rate 
b(z) verify 

0<K(z), JlC(z)dz = l, J K{z)e\' A2 dz < oo, (10) 
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b m <b(z,I) <b M , \V x b(z, I)| < U b(z, I), \b(x,I 1 )-b(x,I 2 )\<L 2 \I 1 -I 2 \, (11) 
where b m , &m> L\ and L2 are constants. Finally for equation ([3]) we replace © and ([7]) by 

min J m ) + b(x, I m )) = 0, max [R(x, I m ) + b(x, I M )] = 0, (12) 



d(R + b) 

- j 2 1 aim — /\4 \ 

where K3 and are positive constants 



\R(x,h) - R(x,I 2 )\ < K 3 \h - I 2 \ and - K 4 < V Q J 1 (x,T) < -K^ 1 < 0, (13) 



In both cases, in the limit we expect n(t, x) = or R(x, I) = 0, where n(t, x) is the weak limit of 
n € (t,x) as e vanishes. Since the latter is possible at only isolated points, we expect n to concentrate 
as Dirac masses. Following earlier works on the similar issue |19|. El (U [28], in order to study n, we 

make a change of variable n e (t, x) = e ^ . It is easier to study the asymptotic behavior of u £ instead 
of n e . In section [5] we study the asymptotic behavior of u t while e vanishes. We show that u e , after 
extraction of a subsequence, converge to a function u that satisfies a constrained Hamilton-Jacobi 
equation in the viscosity sense (see [3J I2TH [TBI |2"2"] for general introduction to the theory of viscosity 
solutions). Our main results are as follows. 

Theorem 1.1. Assume (0)-|^). Let n t be the solution of (Q])-^, and u t = eln(n e ). Then, after 
extraction of a subsequence, u e converges locally uniformly to a function u G C((0, 00) x W 1 ), a viscosity 
solution to the following equation: 

' d t u= |Vu| 2 + R(x,I(t)), 
maxu(t,x)=0, V£ > 0, 



I e (t) — >I(t) a.e., i/}(x)n(t,x)dx = I(t) a.e.. (15) 

e^O J 

In particular, a.e. in t, supp n(t, •) C {u(t, •) =0}. Here the measure n is the weak limit of n t as e 
vanishes. If additionally (u®) e is a sequence of uniformly continuous functions which converges locally 
uniformly to u° then u G C([0,oo) x M rf ) and u(0, x) = u°(x) in M. d . 

Theorem 1.2. Assume [B\)- il3\) . and (u®) e is a sequence of uniformly bounded functions in W 1 ' 00 
which converges locally uniformly to u°. Let n e be the solution of (0)-|^p, and u e = eln(n e ). Then, 
after extraction of a subsequence, u e converges locally uniformly to a function u G C([0, 00) x R d ), a 
viscosity solution to the following equation: 

' d t u = R(x, I(t)) + b(x, I(t)) J K(z)e Vu - z dz, 

< max u(i, x) = 0, Vt > 0, (16) 
^k(0, x) = u°(x), 

I £ (t) — >I(t) a.e., n(t,x)dx = I(t) a.e.. (17) 

e^O J 

In particular, a.e. in t, supp n(t, •) C {u(t, •) = 0}. As above, the measure n is the weak limit of n e 
as e vanishes. 



3 



These theorems improve previous results proved in [191 E 129] in various directions. For the 
case where mutations are described by a Laplace equation, i.e. ([I])-©, Theorem 11.11 generalizes the 
assumptions on the initial data. This generalization derives from regularizing effects of Eikonal hamil- 
tonian (see |26[ [Tl [2]). But our motivation is more in the case of equations ([3])-(j3]) where mutations 
are described by an integral operator. Then we can treat cases where the mutation rate b{x, I) really 
depends on x, which was not available until now. The difficulty here is that Lipschitz bounds on the 
initial data are not propagated on u t and may blow up in finite time (see [HI EJ [15] for regularity 
results for integral hamiltonian) . However, we achieve to control the Lipschitz norm by — u e , that goes 
to infinity as \x\ goes to +oo. 

We do not discuss the uniqueness for equations (fl~4l) and (fl~6|) in this paper. The latter is studied, 
for some particuler cases, in J8[[7]. 

A related, but different, situation arises in reaction-diffusion equations as in combustion (see [U EJ 
HU[ [2"T| [23[ [30] ) . A typical example is the Fisher-KPP equation, where the solution is a progressive 
front. The dynamics of the front is described by a level set of a solution of a Hamilton-Jacobi equation. 

The paper is organized as follows. In section [2] we state some existence results and bounds on n e 
and I e . In section [3] we prove some regularity results for u e corresponding to equations ([I])-©. We 
show that u t are locally uniformly bounded and continuous. In section [4] we prove some analogous 
regularity results for u t corresponding to equations ([3])- ([!]). Finally, in section [5] we describe the 
asymptotic behavior of u e and deduce the constrained Hamilton-Jacobi equation (|14j) -(|15 p . 



2 Preliminary results 

We recall the following existence results for n e and a priori bounds for I t (see also [H 117]). 

Theorem 2.1. With the assumptions $5§-$3$, and I m — Ce 2 < I e (0) < Im + Ce 2 , there is a unique 
solution n e £ C(M + ; L 1 (M d )) to equations (QP-fDJ and it satisfies 

I' m = lm- Ce 2 < I e {t) < I M + Ce 2 = l' M , (18) 

where C is a constant. This solution, n e (t,x), is nonnegative for all t > 0. 

We recall a proof of this theorem in Appendix [A] We have an analogue result for equations ([3])-([4]): 

Theorem 2.2. With the assumptions [W\) - [T3\) . and I m < I e (0) < Im, there is a unique solution 
n e £ COEt^L 1 nL°°(R d )) to equations and it satisfies 

I m < hit) < Im- (19) 

This solution, n e (t,x), is nonnegative for all t>0. 

This theorem can be proved with similar arguments as Theorem 12. 11 A uniform BV bound on I e {t) 
for equations ([I])-© is also proved in [8]: 

Theorem 2.3. With the assumptions we have additionally to the uniform bounds [T8\) . the 

locally uniform BV and sub-Lipschitz bounds 

j t h{t) >-eC + e^ f ^( x )n° e (x)^^dx, (20) 
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j t Qe{t) > -Ct + J {l + ^x))n° e (x)^^-dx, (21) 

where C and L are positive constants and g € (t) = J Rd n e (t,x)dx. Consequently, after extraction of a 
subsequence, I e (t) converges a.e. to a function I{t), as e goes to 0. The limit I(t) is nondecreasing as 
soon as there exists a constant C independent of e such that 

^x)n%x)^l>-Ce^. 

We also have a local BV bound on I e (t) for equations (|3])-([4|): 

Theorem 2.4. With the assumptions ^\)- p3\) . we have additionally to the uniform bounds flgj) . the 
locally uniform BV bound 

| f , (t) >_ C - + e ^| n o (x) ^. J ?) + ^.i?)^ (22) 

r T d 

J \—I e (t)\dt < 2C'T + C", (23) 

where C , C" and V are positive constants. Consequently, after extraction of a subsequence, I t (t) 
converges a.e. to a function I(t), as e goes to 0. 

This theorem is proved in Appendix IB1 

3 Regularity results for equations (p])-(|2]) 

In this section we study the regularity properties of u t = e In n e , where n e is the unique solution of 
equations ([I])-©. We have 

dtn e = -dtu t e~ , Vn t = -Vu e e~, An e = -An e + \Vu e \ )e~ . 

e e e e z ' 

Consequently u e is a smooth solution to the following equation 



J d t u t - eAu t = \Vu t \ 2 + R(x, I e (t)), x £ R, t > 0, 
\u e (t = 0) = elnn°. 

We have the following regularity results for u e . 



(24) 



Theorem 3.1. Let T > be given. Then u € < D 2 , where D = B + (A 2 + K<i)T , and we define 
v e = \/2D 2 — u e . With the assumptions for all to > v e are locally uniformly bounded and 

Lipschitz in [to , T] x M. d , 

\Vv t \<C(T) + ^=, (25) 

where C(T) is a constant depending on T, K\, K2, A and B. Moreover, if we assume that (u®) £ is a 
sequence of uniformly continuous functions, then u t are locally uniformly bounded and continuous in 
[0,oo[xM d . 

We prove Theorem 13.11 in several steps. We first prove an upper bound, then a regularizing effect 
in x, then local L°° bounds, and finally a regularizing effect in t. 
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3.1 An upper bound for u t 

From assumption ([9]) we have u®(x) < — A\x\ + B. We claim that, with C = A 2 + K2, 

u e (t, x) < -A\x\ + B + Ct, Vt > 0. (26) 
Define (f>(t, x) = -A\x\ + B + Ct. We have 

9t0 - eA0 - IV0I 2 - R(x, Ut)) >C + e A ^ d ~ ^ - ,4 2 - if 2 > 0. 

Here is an upper bound for /) according to (|8|). We have also </>(0, 2) = — A\x\ + B > u®(x). 
So (p e is a super-solution to ([2^1 and (f26|) is proved. 



3.2 Regularizing effect in space 

Let u = f(v), where / is chosen later. We have 

d t u = f'(v)d t v, d x u = f'(v)d x v, Au = f'(v)Av + f"(v)\Vv\ 2 . 
So equation (p^|) becomes 



dtv — eAv 



/"(») ^ ,,, , 
'-m +S{v] 
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(27) 



Define p = Vt>. By differentiating (|27p we have 



%>; - eApi - 2 



We multiply the equation by and sum over i: 



. f'(v) 



J'»(y) f"(v) 2 



+ /» 



/v 2 



i>j + — 



\Vv\ 2 Pi 

1 ai? 



Vv ■ V- 



r» /» 2 



First, we compute Yli(^Pi)Pi- 



|p|AH + |V|p|| 2 -^|Vp^ 
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We also have 



ivni 2 = E * E i^i 2 = E i^-i 2 = E i^'i 2 - 



>, j 



It follows that 



E(Aft)ft < |p|A|p|. 



We deduce 



d t \p\ - eA\p\ - 2 



P- V|p| 



f'M l » 



(28) 



Prom (HH) we know that, for < t < T, ?i £ < £>(T) 2 , where Z?(T) = y/B + CT. Then we define 
/(u) = — w 2 + 2D 2 , for u positive, and thus 



D(T) < v, 



f(v) = -2v, and 



1 1 

< 



fi{vy 2v 2D 



f"( v ) = -2, and 



/» | 

/» 2 ' 2w 2 " 2D 2 



1 1 

< 



f'iy) = o, 



+ /» 



From (|28p . assumption ([8]) and these calculations we deduce 



2 + e^ > 2. 



d\p\ 
dt 



eA|p| - 2 



e m +f{v) 



p ■ V\p\ + 2\pf 



2D^ ~ ~23 ~ °" 



Thus for 0(T) large enough we can write 



d\p\ 
dt 



eA\p\ -2 



e 7>) +/(?;) 



p- v|p| + 2(|p| -ef < o. 



(29) 



Define the function 
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y(t,x) = y(t) 



1 



2y/t 



Since y is a solution to (j29j) . and g/(0) 



oo and \p\ being a sub-solution we have 



\p\(t,x) < y(t,x) 



1 



+ 6. 



2y/t 



Thus for v e = \/2D 2 — u € , we have 



Vv e \(t,x) < 



1 



+ 0(T), 0<t<T. 



(30) 



2Vt 



3.3 Regularity in space of u € near t = 

Assume that u® are uniformly continuous. We show that u € are uniformly continuous in space on 



For 5 > we prove that for h small |u e (t,x + h) — u e (t, x)\ < S. To do so define w e (t,x) = 
u e (t,x + h) — u e (t,x). Since n° are uniformly continuous, for /i small enough |u) e (0,x)| < |. Besides 
w e satisfies the following equation: 

d t w e (t,x) - eAw e (t,x) - (Vu e (t,x + h) +Vu e (t,x)) ■ Vw e (t,x) = R(x + h,I e {t)) - R{x,I e (t)), 
Using assumption ([8|) we have 



So for h small enough \u e (t, x + h) — u e (t,x)\ < 5 on [0, T] x Mr. 
3.4 Local bounds for u e 

We show that u € are bounded on compact subsets of ]0, oo[x]R d . We already know from section [3TT1 that 
u e is locally bounded from above. We show that it is also bounded from below on C = [to, T] x B(0, R), 
for all R > and < t < T. 

From section [3TT1 we have u e (t,x) < — A\x\ + B + CT. So for e < eo, eo small enough and R large 
enough 



[0, T] x M d . 



dtw e (t, x) — eAw e (t, x) — (Vu e (t, x + h) + Vu e (t, x)) ■ X7w e (t, x) < K2\h\. 



Therefore by the maximum principle we arrive at 



maxlwJt, x)\ < max\w e (0, x)\ + K2\h\t. 

R d R d 




We have also from (fTgj) that 
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f e u -f dx > Ik, 

JM. d WM 

We deduce that for all < e < eo and R large enough 

/ ^dx > Jjk 

■/|a:|<.R 2 ^M 

Therefore there exists ei > such that, for all e < ei 

Eb G M d ; |zo| < «e(*,«o) > -1, thus u e (* 5 z ) < \/2D 2 + 1. 
From Section 13,21 we know that v e are locally uniformly Lipschitz 

\v e (t,x + h)-v e (t,x)\ < (C(T) + —=)\h\, 
Thus for all (t, x) G C and e < e\ 

v e (t, x) < Efy, T, R) := V / 2D 2 (T) + 1 + 2(C(T) + 1 r^=)R. 

It follows that 

u e (t, x) > 2D 2 (T) - E 2 (t ,T,R). 
We conclude that u e are uniformly bounded from below on C. 

If we assume additionally that u® are uniformly continuous, with similar arguments we can show 
that u e are bounded on compact subsets of [0, co[x]R rf . To prove the latter we use uniform continuity 
of u 6 instead of the Lipschitz bounds of v e . 

3.5 Regularizing effect in time 

From the above uniform bounds and continuity results we can also deduce uniform continuity in time 
i.e. for all rj > 0, there exists 9 > such that for all (t,s,x) G [0,T] x [0,T] x B(0, f), such that 
< t — s < 9, and for all e < eo we have: 

\u e (t, x) — U e (s, x)\ < 7]' . 

We prove this with the same method as of Lemma 9.1 in [3] (see also for another proof of this 
claim). We prove that for any r/ > 0, we can find positive constants A, B large enough such that, for 
any x G B(0, ^) and for every e < eo 

u e (t, y) - u e (s, x) < r] + A\x - y\ 2 + B(t - s), for every (t, y) G [0, T] x B(0, R), (31) 

and 

u € {t, y) - u € {s, x)>-rj- A\x - y\ 2 - B(t - s), for every (t, y) G [0, T] x B(0, R). (32) 
We prove inequality (f3T|) . the proof of ([32]) is analogous. We fix (s,x) in [0,T[xB(0, |). Define 

£(i, y) = u e {s, x) + rj + A\y - x\ 2 + B(t - s), (t, y) G [s, T[xB(0, R), 
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where A and B are constants to be determined. We prove that, for A and B large enough, £ is a 
super-solution to J23D on [s, T] xB(0, i?) and£(i,y) > u e {t,y) for (t,y) G {s} xB(0, R)U[s, T] xdB(0, R). 

According to section 13.41 u t are locally uniformly bounded, so we can take A a constant such that 
for all e < eo, 

. 8 II u e llL°°([0,T]xB(0,fl)) 

" R? ' 

With this choice, £{t,y) > u e (t,y) on dB(0,R) x [0,T], for all 77, B and x G B(0, f ). Next we prove 
that, for A large enough, £(s,y) > u e (s,y) for all y £ B(0,i?). We argue by contradiction. Assume 
that there exists rj > such that for all constants A there exists yA,t E B(0,i?) such that 



it e (s, - « e (s, x) > r? + A|y^ ie - x\ 2 . (33) 

It follows that 

[2M 

\yA,e-x\ < y — , 

where M is a uniform upper bound for || u e ||l°°([o,t]xB(o,.r))- Now let A — > 00. Then for all e, 
|y^4 )£ — x| — > 0. According to Section 13.31 u t are uniformly continuous on space. Thus there exists 
h > such that if |y^ je — x\ < h then \u e (s, yA,e) — u € (s,x)\ < |, for all e. This is in contradiction with 
(|33l) . Therefore £(s,y) > u € (s,y) for all y G B(0, i2). Finally, noting that i? is bounded we deduce 
that for £> large enough, £ is a super-solution to (|24|) in [s,T] x B(0,J2). Since u e is a solution of (f24|) 
we have 

«e(*,2/) < f(t,y) = u £ (s,a;) + r? + A\y - x\ 2 + B(t- s) for all (t,y) G [s,T] x B(0,i2). 

Thus ()3ip is satisfied for t > s. We can prove (j32[) for t > s analogously. Then we put x = y and 
we conclude. 



4 Regularity results for equations (S])-([4]) 

In this section we study the regularity properties of u t = e In n e , where n t is the unique solution of 
equations ©-(111) as given in Theorem I2.2L From equation ([3]) we deduce that u e is a solution to the 
following equation 

u e (t,x-\-ez) — u e {t,x) 

d t u e = R(x,I e (t)) + / K(z)b(x + ez,I e )e i dz, x E R, t > 0, . . 

u e (t = 0) = elnn°. 

We have the following regularity results for u e . 

Theorem 4.1. Define v e = n\ = exp(u e ) ; where n e is the solution to equations ^)-§4ty. With the 
assumptions ^)- lll3\) . and if we assume that (u®) 6 is a sequence of uniformly bounded functions in 
W l,oc , then u e are locally uniformly bounded and Lipschitz in [0, oo[xR d . 

As in section [3] we prove Theorem 14.11 in several steps. We first prove an upper and a lower bound 
on u e , then local Lipschitz bounds in space and finally a regularity result in time. 
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4.1 Upper and lower bounds on u e 

From assumption ([9]) we have u®(x) < —A\x\ + B. As in section [3TT1 we claim that 

u e (t,x) < -A\x\ +B + Ct, Vt > 0. (35) 
Define v(t,x) = -A\x\ +B + Ct, where C = b M J K(z)e A ^dz + K 2 . Using © and (HH) we have 



A \*\dz > 0. 



d t v - R(x, I e (t)) - J K(z)b(x + ez, I e )e" ( *'' B+ " ) "^tfe >C-K 2 -b M j K(z)e 

We also have u(0, x) = — -A|x| + -B > u°(x). So u is a supersolution to (f3~i|) . Since ([3]) verifies the 
comparison property, equation (|34|) verifies also the comparison property, i.e. if u and v are respec- 
tively super and subsolutions of (f34"|) then v < u. Thus (f35j) is proved. 



To prove a lower bound on u e we assume that u® are locally uniformly bounded. Then from equation 
(j34"|) and assumption © we deduce 

d t u e (t,x) > -K 2 , 

and thus 

u e (t,x) > -\\u° e \\ L oo(B(o,R)) ~ K 2 t, Mx G B(0,i?). 
Moreover, \Vu®\ being bounded, we can give a lower bound in R. d 

uJt,x) > infit°(0) - || Vn°|| L oc|x| - K 2 t, Vx G M. d . (36) 
4.2 Lipschitz bounds 

Here we assume that u e is differentiable in x (See [E]). See also Appendix [C] for a proof without any 
regularity assumptions on u e . 

Let p e = Vit e ■ x, where x is a fixed unit vector. By differentiating (f34"|) with respect to x we obtain 

it e (t,x-\-ez)—u e (t,x) 

K{z)Vb(x + ez, I e ) ■ x e = 
+ / g(^)6(x + ez, J e ) Pe( *' X + €Z) - P ^ X ) e ^ + T^ dz. 



Thus, using assumptions dH) and (jlip . we have 



/ite(t,a; + ez) — u^(t,x) 
K(z)b(x + ez,I e )e i (37) 

+ / g(*)6(s + ez, J e ) Pe( *' X + ez) ~ P ^ X) dz . 
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Define w e (t,x) = p e (t,x) + L lUe (t,x) and A € (t,x,z) = u ^ x+ez ) From ^ and ^ we 

deduce 

d tWe - K 2 (l + Ll )-j K(z)b(x + ez, l£) Mt^ + ez)-w e (t,x) eAeit>x>z)dz 

< 2Li J K(z)b{x + ez,I € )e A ^ x > z) dz 

= Li y tf(z)&(a; + ez, I e )e A ^ x ^(2 - A e (t,x,z))dz 

< Lib M e, 

noticing that e is the maximum of the function g(t) = e (2 — t) in M. Therefore by the maximum 
principle, with C\ = K 2 (l + L\) + Lib^e, we have 

w; e (i, x) < Cii + max w € (0, x). 

It follows that 

Pe(t,x) < dt+ || Vu° || L0O + J L 1 ( J B + Ct)+ J L 1 (||Vu°|| i o O |x|+ J Fs:2i-«°(a; = 0)) (38) 
= C 2 t + C 3 \x\ +C 4 , 

where C 2 , C3 and C4 are constants. Since this bound is true for any \x\ = 1, we obtain a local bound 
on |Vu e |. 

4.3 Regularity in time 

In section I4T21 we proved that u 6 is locally uniformly Lipschitz in space. From this we can deduce that 
dtu e is also locally uniformly bounded. 

Let C = [0, T] x B(xo, R) and Si be a constant such that || u e \\l oc (C)< Si f° r an e > 0. Assume that 
R! is a constant large enough such that we have u e (t,x) < —Si in [0, T] x M. d \B(xo, R'). According to 
([35]) there exists such constant R'. We choose a constant S 2 such that || Vu e ||L oo ([o,T]xB(a;o,-R')) < 
for all e > 0. We deduce 

K{z)b{x + ez, I t )e i 1 

\x+ez\<R' 

+ 1 \x+Ez\>R')dz 

<K 2 + b M J K{z)e s ^ z h\ x+ez \ <Rl dz + b M j K(z)l\ t +ez\>R'dz 
<K 2 + b M (l + j K(z)e s '^dz). 

This completes the proof of Theorem 14.11 
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5 Asymptotic behavior of u e 



Using the regularity results in sections [3] and HJ we can now describe the asymptotic behavior of u e 
and prove Theorems 11.11 and 1 1 . 2 1 Here we prove Theorem ll.il The proof of Theorem 1 1.2 1 is analogous, 
except the limit of the integral term in equation (|16p . The latter has been studied in |X9|, [T2| [JJ [29] . 

Proof of theorem \1. 11 step 1 (Limit) According to section [31 u 6 are locally uniformly bounded and 
continuous. So by Arzela-Ascoli Theorem after extraction of a subsequence, u e converges locally uni- 
formly to a continuous function u. 

step 2 (Initial condition) We proved that if u® are uniformly continuous then u e will be locally 
uniformly bounded and continuous in [0, T] x Thus we can apply Arzela-Ascoli near t = as well. 
Therefore we have u(0,x) = lim u e (0,x) = u (x). 

e— >0 

step 3 (maxii = 0) Assume that for some t,x we have < a < u(t,x). Since u is continuous 

u(t,y) > | on B(x,r), for some r > 0. Thus we have n e (t,y) — > oo, while e — > 0. Therefore I e {t) — > oo 
while e — > 0. This is a contradiction with (1181). 



To prove that max u(£, x) = 0, it suffices to show that lim n e (t, x) ^ 0, for some x 6 R d . From ([26 



we have 

u e (t,x) < -A\x\ + B + Ct. 

It follows that for M large enough 



lim 

From this and (fT8l) we deduce 



f f -A\x\+B+Ct 

i I n e (t,x)dx < lim / e J = 0. (39) 

°ib|>M <E - >0 J|x|>M 



lim / n e (t, x)dx > 

e ^°J\x\<M 



4>M 

If u(t, x) < for all \x\ < M then lim e J = and thus lim I, ,^ tf n e (t,x)dx = 0. This is a 
contradiction with ([39]) . It follows that max u(t,x) = 0, Vt > 0. 



step 4 (supp n(t, ■) C •) = 0}) Assume that u(to, xq) = —a < 0. Since u e are uniformly contin- 
uous in a small neighborhood of (to, xo), (t,x) £ [t^ — 5, to+5] x~B(xq,5), we have u e (t,x) < — § < for e 

small. We deduce that w . x . lX] „-n, x ^n(t,x)dtdx = ,, w , x , , r,„i,/ ^ lim e J dtdx = 0. 

Therefore we have supp n(i, •) C {«(£, •) =0}. 

step 5 (Limit equation) Finally we recall, following [8], how to pass to the limit in the equation. 
Since u € is a solution to ([24"[) . it follows that <p e (t,x) = u e (t,x) — f* R(x,I e (s))ds is a solution to the 
following equation 

d t <f) t {t,x) - eA(j) e (t,x) - |V(/> e (t,x)| 2 - 2V(/) € (t,x) • f VR(x, I € (s))ds 

Jo 

r-t ft 



e [ AR(x,I e (s))ds + | f V R(x, I £ (s))ds\ 
Jo Jo 
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Note that we have I e (s) — > I(s) for all s > as e goes to 0, and on the other hand, the function 
R(x, I) is smooth. It follows that we have the locally uniform limits 



lim f R(x,I e (s))ds = [ R(x,I(s))ds, 
€ ^°Jo Jo 

lim/ VR(x,I e (s))ds = [ VR(x, I(s))ds, 
e ^°Jo Jo 

lim I AR(x,I e (s))ds = [ AR(x, I(s))ds, 
e ^°Jo Jo 



for all t > 0. Moreover the functions f Q R(x,I(s))ds, f VR(x,I(s))ds and J* AR(x,I(s))ds are 
continuous. According to step 1, u t (t,x) converge locally uniformly to the continuous function u(t,x) 
as e vanishes. Therefore <p t (t,x) converge locally uniformly to the continuous function (j)(t,x) = 
u(t,x) — Jq R(x, I(s))ds as e vanishes. It follows that (f>{t,x) is a viscosity solution to the equation 

d t (f>(t,x) - \Vc/)(t,x)\ 2 -2V<f>(t,x) • I VR(x,I{s))ds 

Jo 

= | f VR(x,I)d. 
Jo 



si 2 . 



In other words u(t, x) is a viscosity solution to the following equation 

d t u(t, x) = \Vu(t, x)\ 2 + R(x, I(t)). 

A Proof of theorem 12.11 

A.l Existence 

Let T > be given and A be the following closed subset: 

A = {u e C([0,T],L 1 (M rf )), u > 0, || u(t,-) \\ L i< a}, 

where a = [J nP e dx) e~ ?~ . Let $ be the following application: 

$ : A -» A 

where u is the solution to the following equation 



d t v - eAv = ^R{x, J«(t)), i 6 I, i > 
u(< = 0) = n° e . 



□ 



(40) 

/„(t) = / if>(x)u(t,x)dx, (41) 
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and R is defined as below 

'R(xJ) if !f<K2I M , 



R(x,I) 



R(x,2I M ) if 2I M <I, 
R{x,if) if I< l -f. 



We prove that 

(a) $ defines a mapping of A into itself, 

(b) $ is a contraction for T small. 

With these properties, we can apply the Banach-Picard fixed point theorem and iterate the con- 
struction with T fixed. 

Assume that u G A. In order to prove (a) we show that v, the solution to (]40p . belongs to A. By 
the maximum principle we know that v > 0. To prove the L 1 bound we integrate (|40p 



— / vdx = [ —R(x, I u {t))dx < - max R(x, I u (t)) [ vdx < — - / vdx, 
dt J J e e X £R d J e J 



and we conclude from the Gronwall Lemma that 



v II T.i < { I n®dx ) e ? = a. 



Thus (a) is proved. It remains to prove (b). Let u±, u 2 G A, v\ = §{u\) and V2 = &(u 2 ). We have 

d t (vi - v 2 ) - eA(t>i - v 2 ) = - [(vi - v 2 )R(x,I Ul ) + v 2 (R(x,I Ul ) - R(x,I U2 ))] . 

Noting that || v 2 \\ L i< a, and \R(x,I Ul ) - R(x,I U2 )\ < Ki\I Ul - I U2 \ < Kiip M \\ ui - U2 Wl 1 we 
obtain 

d ,. ,. . K 2 1aK\^M ,, ,, 

37 II v l ~ v 2 ZA< \\Vl- V 2 \\ L l H || Ui - U 2 \\ L i . 

dt € € 



Using v\(0, •) = v 2 (0, •) we deduce 



II n - v 2 \\l^li< — y 2 — ( e £ ~ 1 ) II Ul ~ M2 ■ 

K 2 T K 2 T K 

Thus, for T small enough such that e e (e e — 1) < iKl ^^ j n o > ^ is a contraction. Therefore $ 
has a fixed point and there exists n e G A a solution to the following equation 

(d t n e - eAn e = ?fR(x, I{t)), x G R, < t < T, 
\n e (t = 0) = n°. 



I(t) = / ip(x)n e (t,x)dx, 

With the same arguments as IA.2I we prove that kf- < I(t) < 21m and thus n e is a solution to 
equations ©-© for t G [0,T]. We fix T small enough such that e^(e^ - 1) < ^k^mIm ' Then 
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we can iterate in time and find a global solution to equations ([I])-©. 



Applying the maximum principle to the equation we deduce that n e is nonnegative. 

A. 2 Uniform bounds on I e (t) 
We have 



dl d f f 1 

— — = — / ijj(x)n e (t,x)dx = e / ■ip(x)An e (t, x)dx H — / ip{x)n e {t,x)R{x,I e {t))dx. 
dt at .had ./rod e 



We define = XL • tfi £ WJf^QR ), where xl is a smooth function with a compact support such 
that xl\b{o,l) = 1, Xl|r\b(o,2L) = 0. Then by integration by parts we find 



i[>l(x) An t (t, x)dx = I AijjL(x)n e (t, x)dx. 
As L — > oo, ipL converges to ip in W 2 '°°(R d ). Therefore we obtain 

lim / Ai^L(x)n e dx = / A^(x)n e dx, 



L— >oo 



lim 

L— >oo 



/ ipL{x)An e (t,x)dx = / ip(x)An € (t,x)dx. 



From these calculations we conclude 



It follows that 



di f i 

— - = e / Aip(x)n e (t, x)dx H — / ip(x)n e (t,x)R(x, I £ (t))dx. 
dt ./rod e 



C 1 <i/ C 1 

e- — I e H — I e min R(x,I e ) < < e I e H — I e max R(x, I e ). 



Let C = . As soon as I e overpasses 1m + Ce 2 , we have R(x,I e ) < = — anc ^ thus 



die 

becomes negative. Similarly, as soon as I e becomes less than I m — Ce 2 , ^ becomes positive. Thus 
(fT8|) is proved. 



B A locally uniform BV bound on I e for equations ([3])- ([4 

In this appendix we prove Theorem 12.41 We first integrate (|3|) over M. d to obtain 

d I t (t) = - f n e (t,x)(R(x,I € (t)) + b(x,I e (t)))dx. 



dt cw e . 

Define J e (i) = ^I t {t). We differentiate J e and we obtain 
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j t J e (t) = ~ J e (t) J n e {t,x) d ^ R ^ h \ x,I e (t))dx 

+ {R(x,Ie) + b(x,I e ))[n e (t,x)R(x,I e ) + J K e (y - x)b(y,I e )n e (t,y)dy]dx. 

We rewrite this equality in the following form 

j t J e {t) = ^J e (t) J n £ (t, x) 8 ( R d + 6) (x, I e (t))dx + ^ J n e (t,x)(R(x,I)+b(x,I)) 2 dx 
+ ^ y J K € (y-x)(R(x,I)-R(y,I))b(y,I e )n e {t,y)dydx 
+ ~2 f f K z(y ~ x )(b( x ,I) -b(y,I))b(y,I e )n e (t,y)dydx. 

It follows that 

j t Je(t)>^Je(t) J n £ (t, x) ™ R + h > (x, I e (t))dx + ^ y n e (t,x)(i?(a:,/) + 6(x,/)) 2 dx 

+ 62, I € )n € (t, x + ez)dzdx 



K 2 + bM L 



>^J e (t) J n £ (t,x)-^^-(x,I e (t))dx + ^ j n e (t,x)(R{x,I) + b{x,I)) 2 dx 



where C\ is a positive constant. Consequently, using (|T3j) we obtain 



^(*))-<^-^W*))- 

ar e e 



From this inequality we deduce 



(J.(t))-<^ + (Je(0))_e-^. 
With similar arguments we obtain 

(Je(t)) + >-^ + (Je(0)) + e-^. 

Thus (|22p is proved. Finally, we deduce the locally uniform BV bound ([23 



<I M -I m + 2C'T + 0{1). 
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C Lipschitz bounds for equations ([3])- ([4]) 

Here we prove that u t are locally uniformly Lipschitz without assuming that the latter are differn- 
tiable. The proof follows the same ideas as in section I4.2L 

Let c = From ([Ml) we have 

d t (u e (t,x + h) - u e (t,x) + ch (2u e (t, x + h) - u e (t,x)) - (1 + 2ch)R(x + h,I e ) + (1 + ch)R(x,I e ) 

/lie (t,x J rh-\-ez) — u £ (i ,x-\-h) j Ug (t ,x + ez) — ii e (t,x) 

K(z)b(x + h + ez, I e )e J — / K(z)b(x + ez, I e )e ' dz 

, j u £ (t,x J rh-\-ez) — u^(t,x-\-h) j it e (t,x-\-ez) — u £ {t,x) . 

+ c/i( / K(z)2b(x + h + ez,I e )e J - / K(z)b(x + ez, I € )e i dz) 

Denne Q = ^ = ? ^ x) = x + h) _ u ^ x) and ^ x) = 

u e (t,x+h)-u c (t,x) _j_ ^ Using the convexity inequality 

e < e « + e /3( ( g_ a ) ; 

we deduce 

hd t w e (t, x) - (1 + 2ch)R(x + /i, I e ) + (1 + ch)R{x, I e ) 
< / K(z) K * + * + «, + ^ - „))„, - / K(m , + «, 

+ c/i( y 2K(z)b(x + h + ez,I t )e /3 dz- J K(z)b(x + ez, I t )e a dz) 
<j K{ z m x + h + ez, Q -b(x + ez, Q )e«dz 
+ y + h + ez, I e )eP (/? - a + c/ ^ ( *' - + ^ ~ A ^ x) )dz 

+ chj K(z)b(x + h + ez,I e )e t3 (2-2f3 + a)dz-ch J K(z)b(x + ez, I e )e a dz. 
From assumptions ([8]) and (jlip it follows that 

d t w e (t,x)< / K(z)6(x + /i + ez,/ £ )e /3 -^ ^ ^-^dz 

+ + 3cK 2 + y K(z)(cb M e P (2 - 2(3 + a) + (L x b M -cb m )e a )dz. 

Notice that 

o&M e^(2 - 2/3 + q) + (Li&m - cb m )e a = cb M e P {2 - 2/3 + a) - Li& M e a , 

is bounded from above. Indeed if we first maximize the latter with respect to (3 and then with respect 
to a we obtain 
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cb M e f3 {2 -2(3 + a)- L x b M e a < 2cb M e% - L x b M e a < 

We deduce 

/ \ [ tj-, m, , T x f,W e (t,X + €Z) — W e (t,x) , _ 

d t w e (t,x)< / K(z)b{x + h + ez,I e )e f) -^ - -^—t-dz + G, 

where G is a constant. Therefore by the maximum principle, ([35 1) and (|36p we have 

w e (t, x) < Gt+ || Vu° e \\ L oo -2cA\x + h\ + 2cB - cu° t (x = 0) + c || Vu° \\ L °° \x\ 
Using again ([35]) and ([36]) we conclude that 



u ^ x + h l u ^ x ~) <{G + 2cK 2 )t + c{-A+ || Vu° \\loc ) (|x| + 2\x + h\) (42) 



+ 3cB+ || Vu° e || L oo -3cinf = 0). 
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